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ABSTRACT 

A method of digital computer oriented analysis of a 
Mason flow graph is presented. Using this method in connec- 
tion with a computer program implementing the algebraic 
manipulation language FORMAC can be a tool in expanding large 
determinants with symbolic polynomial elements and thereby 
finding the characteristic equation and transfer functions of 
a linear time-invariant system in symbolic form. 

Also a new approach to the root locus method is shown, 
using a FORMAC program. The advantages over the conventional 
root locus method are dciseueeaae 

Areas of possible future use of FORMAC in algebraic 


problems of control theory are discussed. 
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fen eRO DUCTION 


ye THE NATURE OF THE PROBLEM 

Consider a linear time-invariant system whose parameters 
are known only symbolically. In the analysis of the system 
often it becomes necessary to derive certain system functions 
which are polynomials in the Laplace variable s. 

The "manual" derivation of these functions, even for 
small systems, becomes a task which is not only time consum- 
ing but unattractive, and although conceptually simple, is 
also a source of errors which are difficult to detect. 

The enormous work involved in the deriving of these 
functions for medium or large systems can be the reason for 
not attempting a solution to the problem in general algebraic 
fOr . 

There are also possible other algebraic problems where 
the "manual" derivation is tedious. 

The problem, therefore, is to devise methods of using 


the digital computer to perform those algebraic manipulations. 


Bi. METHOD OF SOLUTION 
The first specific problem to be considered will be the 
evaluation of the characteristic equation of a linear system. 
To solve this problem it will be necessary to evaluate 
the determinant of a square matrix whose elements are generally 
polynomials in the variable s, derived from the differential 
equations of the system by Laplace transform. 
Two different methods of evaluating a determinant are 


applied here. The first method is based on the theory of 
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graphs, specifically on the properties of a Mason flow graph 
(1). It will be assumed that the reader is familiar with the 
Mason gain formula. 

Other graph methods are those of Coates and the modified 
method of Chan and Bapna [Ref. 2]. These methods are not 
based on a agonmel ou graph. A flow graph analysis using a 
Mason flow graph and based on a topological development is 
given by Dunn and Chan [Ref. 3]. 

The flow graph analysis used here is not pure topological, 
because the loops and paths are not found from the topological 
properties of the graph but by search methods. 

A comparison of flow graph techniques and another 
different method by Chan and Mai are found in [Ref. 2]. 

The basic feature of the mechanized graph analysis used 
here is the process of identifying closed unidirectional paths 
called “loops" among the connections of the graph. Then the 
set of non-touching loops is found which is then used to 
evaluate the determinant. 

The second method used here for the evaluation of a 
determinant is using a fraction-free algorithm. Lipson [Ref. 
4} uses this algorithm in the evaluation of a set of linear 
algebraic equations. 

The basic problem in the evaluation of the determinant 
Det (A) is the identification of those elements in the square 
matrix A whose products form a term in the determinant. "imere 


are n! terms in the determinant of a Square matrix of order 


n; most of them are often equal to zero in a linear system of 
the kind considered in control theory. 

The computer program DTERM (Program 1) shows in its 
algebraic part how products of polynomials can be processed 
without using excessive space (which was the main limitation 
imposed). 

The first, non-algebraic part of DTERM uses the basic 
definition of a determinant to find the non-zero products in 
the determinant by forming all possible permutations of the 
column indices of the elements of the matrix A. This method 
was used only to simulate the presence of a result of the 
flow graph analysis, which was not programmed. 

For matrices of order 5 it was found that the execution 
time for investigation of all 120 possible products of 5 
elements was 6 sec, for order 8 the execution time was 27 sec 
and for order 9 the execution time was 3 min 43 sec. These 
times include also other statements. 

The fraction-free algorithm is given here as a different 
possible method to evaluate the determinant of a matrix, using 
only single letters as symbolic elements followed by a gradual 
replacement by the longer polynomials. The latter process is 
the problematic one because of the enormous space requirement 
which can be generated. It would probably be similar to the 
method used in the algebraic part of DTERM. 

For the other algebraic problems mentioned here short 


example programs in FORMAC, an algebraic manipulation language, 


are given. These programs show the possible use of FORMAC in 


the solution of those problems; they are not optimal. 


Ge THE ALGEBRAIC MANIPULATION LANGUAGE "FORMAC" 

The algebraic manipulation must be performed with the 
aid of a special tool: an algebraic manipulation language. 
Here "FORMAC" was used for this purpose. 

FORMAC stands for FORmula MAnipulation Compiler. It was 
developed by the IBM Corporation as "an experimental system to 
provide engineers with the capability of handling formal 
mathematical expressions on a computer"([5]. 

Examples of application of FORMAC are given in'®) ana (7) 
Reference to bibliographies is made in (8) | 

The most recent form of FORMAC is available in the PLl- 
Version. The language PL/1 is a proper subset of FORMAC, it 
is therefore possible to use the entire PL/1 capability toget- 
her with FORMAC. The basic reference manual for PL/1 (9) 


(6) at this 


together with the respective manual for FORMAC 
time (1970) is all the user must be referred to. The average 
FORTRAN user can however in a few days acquire effective know- 
ledge of both PL/1 and FORMAC in order to communicate with the 
FORMAC language. 

FORMAC enables the user to analyze and synthesize general 
algebraic expressions. 

The work with FORMAC proceeds basically in such a way 
that the user writes a program in the FORMAC language, follow- 


ing the route of analysis of normal “manual. “algebraic 


manipulation. 
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The capabilities of FORMAC can be divided into the 


following groups: 


- analysis: 


-synthesis: 


test for mathematical identity, finding 
the coefficient of an expression inside 
a parent expression, finding the high- 
est or lowest power of one expression 
inside a parent expression, finding 

the numerator or denominator of a 
rational expression, finding the main 
operator of an expression (+,-,:,X, 

and so on), finding the number of terms 
of an expression (factors or summands), 
selection of a specified term of an 
expression. 

assignment of an algebraic expression 
to a FORMAC variable, user defined 
function processing, conversion routines, 
well developed transit between "normal" 
values and PL/1 variables and FORMAC 
variables, multiplication, division, 
addition and subtraction, differentia- 
tion with respect to specified and 
unspecified arguments, use of trigono- 
metric and natural functions with 
numeric and symbolic arguments, combin- 


atorial and factorial integer functions. 


Jed 


-support: Simple printout, editing, control over 


level of expansion and evaluation. 
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If. METHOD OF COMPUTER ANALYSIS 
FOR A MASON FLOW GRAPH 


In this chapter a procedure will be described which can 


be implemented using a digital computer. The basic operations 


performed by the computer are list searches and logic oper- 


ations on binary matrices. 


The method produces: 


lists of node numbers which form a specific loop. 
lists of associated branch gains, so loop gains 
can be evaluated. 

a list containing the names of touching loops, so 
sets of touching loops can be identified. 

lists of node numbers which form a path from a 
specified node to another specified node. All the 
lists identify all the paths present between these 
nodes. 

lists of associated branch gains, so the different 
path gains can be evaluated. 

a list of branch point nodes which identify the 


points where the different paths branch. 


First the terms used will be defined, then the method will 


be described, using examples. The computer implementation will 


be discussed. 


A. Bee ENILIONS 


1. The Mason Flow Graph 


Here the definitions and notation of Ward and Strum[1l] 


will be used. 
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The Mason flow graph consists of a set of nodes, 
which are here arbitrarily numbered, and a set of branches, 
which have a direction and a weight associated with them. The 
branches connect the nodes in such a way that at least every 
node is connected with one other node; self loops are not 
present. 

2. Input- and Output-Nodes, Branchpoint Node 

An input node has no incoming branches. Here an out- 
put node has one incoming branch. A branchpoint node is a node 
with at least two outgoing branches. From a branchpoint node 
there are at least two paths leading to a common end node. 

End- and starting nodes are the end- and the Staremnag 
point respectively of paths. 

5. Sihpue fast mca 

The input data are prepared in the form of an input 
list, which is called "LisTi". Each item in LISTl has three 
parts: two node numbers and a branch designator. 


These definitions are illustrated in Figw1 andi 
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EPEC ue 
Nodes 3, 2, 6 are input nodes. This graph contains no 


branchpoint nodes. The other nodes are normal nodes, 


except node 9, which is output node. 





op he ee 


In this graph node 4 is (also) a branchpoint node. 


ile) 


Ese 


Column 1 Column 2 Column 3 
node number k branch designator A node number 1 
from via to 


elements in Col. 1 are referred to as LIST (1,k) 
elements in Col. 3 are rererred t6 as@miSsTL a3 oko 

LIST1l constitutes a complete listing of the graph. 
Columns 1 and 3 contain numbers, Column 2 contains an alpha- 
numeric designator in the form “letter-number", "letter" or 
"letter-letter". 

The order of sabe listing is arbitrary: it is recommend- 
ed to bring the list in at least approximate numerical order in 
Column 1 for improved error detection: 

4. LIST2 and LIST3 

LIST2 1s identical to LIST] except Columns 1 and 3 
are exchanged and then Column 1 in LIST2 is arranged in 
humerical order. LIST3 has the same structure as LISTl, but 
the contents have changed. 

2. Topological Matrices 
a. Loop Matrix 
A loop matrix L is defined. The elements Of GE 
are aj). 
ain = ii for node k being in loop j 
= 0 O.W. 
J=1,2,3,...nl, nl = number of loops 


k=1,2,3,...nn, nn = number of nodes 
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moo O 0 0 0 
Coe bet 0° 0 0 
oo 1 ] ini 


oeoy oO: Cir 


|'d 


Baliga 3 
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b. Path Matrix Pp 


A path matrix P is defined. The elements of P 


are bit. 
ae =i for node k being in path m 
= 0 O.W 
m= 1,273) aaenn, nn = number of nodes 


c. Touching Loops Matrix LT 
A matrix is defined which contains information 
about loops which touch in the Mason flow graph. The elements 
of this matrix LT are Gey: 
ore = 1 for loop k touches loop 7 3"tha eae loops 


k and j have at least one node in 


common. 


k, j = 1,273, en! , gn = number of loops 
d. "Path - Loop Touching” Matrix Pr 
A matrix PL is defined which indicates which 
paths touch which loops. The elements of PL are e,_ 
ees i for loop k touches path m 
= 0 O.W. 


K=1 42 Sse er number of loops 


I 


m=1,2,3,.@eip, np number of paths 
e,. Control Vectors 
Certain control vectors are defined to check on 


Specified conditions. 
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(1) Loop Search Control Vector LC 
The element of LC are Cig: 
= 1 for every input and output node 
= 1 for every other node used at least once 
in the loop search 
= 0 OD. W. 
k=1,2,3,...nn, nn = number of nodes 


(2) eMultiple Node Control Vector LIST1 (4) 
mor (4,k) 


O for node LIST1(1,k) being a single 
node or the mth node in a series of 
m equal node numbers 
= 1 for the first m - 1 node numbers in 
a sequence of m equal node numbers 
k=1,2,3,..-nn, nn = number of nodes 
(3) Branchpoint Control Vector BC 
This vector contains information about nodes 
which are branchpoint nodes. The information is in the form 
of an index k, identifying the branchpoint node in LIST2(1,k). 
BE) e=-.0 for node i being not a branchpoint node 
= ji for node 1 (= LIST2(1,k) being a branch- 
point node. 
6. Logical Operations On Matrices 
a. The Operation @ 
By the ® - product of two matrices A @ B in that 
order of the m x p matrix A = (ai4) and the p x n matrix B = 


(by 4) 1S meant the m x n matrix C = (ci 5) where 
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C.. = asa & b,. ®a & b Oe ke O a. & b 


ij 1) eZ 


i 
: stg & Daa)s a=1,2,...m; j=l1,2,...n 


where & is the logical operation AND with truth table: 


a b a&b 
0 0 0 
0 il 0 
il 0 0 
al 1 1 


and ® is the logical operation OR with the truth table: 


a b c®b 
0 0 0 
il i 1 
0 1 1 
1 0 di 
lore THE LOOP SEARCH 
1. Initialization 


Consider LIST1. Define two temporary control vectors 
Tl and T2 such that for every node there is one binary position 
reserved. Then set for every node k present in LIST1(1l) the 
element T1(k) to 1 and for every node m present in LIST1(3) 
the element T2(m) to m. 

Tl({k)eeeol fore. bIST!] (1 |i AD 


f24{m) = 1. for List (37m 0 
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Now the control vector LC is initially set to 


EeGs="Tl-e T2 


Example: 

ee Sayles ee) es 
3 Bae 
Ai Crud 
4 da 5 
Sae 6 

graph: 





Fig. 4 
Die = “fie 1 leek 26) 
= oe COL 
cea A Oe lL 6 0in OO 1] 


The next step in the initialization is the marking of 
the multiple nodes. Use is made of the control vector LIST1(4). 
The control vector LIST1(4) is set to its value 


according to its definition. (See program FINDLOOP) 
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2. Search For Loops 
a. Finding Successive Nodes 

The Column 1 of LIST1 is called LIST1(1), the 
Column 2 is called LIST1(2). An element in a’ column of LIST1 
will be denoted by LIST1 (column number, element index). 

The search starts in LIST1(1), using LISTER 
as the first search key. The result of this search is an 
index k at which: 

LIST] (3,1) = 4120S e Gas 
Now the search key is changed: 
search key = LIST1(3,k) 
and the search continues, at the top of LIST(1) or at the 
next location. (For short 11ists it may be reasonable to start 
at the top again.) 

When in the continuing search: 

LISTING kK) = LIST a) 
the search key is again changed to: 

search key = LIST1(3,1) 
and the search continues. 

In a graph containing loops, LISTI contaume 
multiple equal node numbers in both Column 1 and 2. 

The question arises how to proceed with the 
search when these multiple node numbers are encountered. The 
control vector LIST1(4) contains information about those nodes 
which have leaving branches in more than one direction. It is 


necessary to follow all possible directed paths through the 


Jags 


graph. LIST1(4) enables the search routine to keep track of 
all nodes where to try to follow also other possible paths. 

If during the search an element LIST1(1,k) is 
encountered whose related control element LIST1(4,k) = 1, then 
the element LIST1(3,k) is selected as the next search key (as 
explained above), but the control vector element LIST1 (4,k) 
is set to "O" before the search continues. 

b. '*Detection Of a Loop 

The element LIST1(1,1) and all the elements 
LISTL(1,n) which are "found" during the search are "marked" in 
a temporary control vector (the max. length is equal to the 
number of nodes). Every time a new search key is selected, a 
eheck 1s made against the control vector. If for the 
Speerere’ node a = LISTT(3,i), a "1" is found in the control 
vector, a loop is detected because a closed directed path had 
been followed, the path being closed at LIST1(3,i). 

c. Recording Of the Loop 

To record the actual sequence of nodes in a loop, 
which is an orderly way of noting the loop nodes, it is 
necessary to store ina list all successively "found" elements 
fiottwiyn), Staremng with the ianitcial valme LISTYY1,ly. This 
is done during the search and the list is abandoned if the 
search is terminated before a loop could be found. But most 
likely the actual sequence is of no consequence. After the 
z#h loop detection the "loop matrix" L is updated by filling 


the j=" row of L with the appropriate entries. 


rs, 


Also for each node touched during the search an 
entry.~as made in.control .vector LC; 

d. Completion Of the Loop Search 

The loop search is basically terminated if by 
the process described in b. above no loop is detected. How- 
ever, before it can be assumed that all loops of the graph 
have been found, the control vector LC must be checked. If 
LC = 1 (contains only "1" ‘'s), all nodes in the graph have 
been checked for belonging to a loop. Otherwise the search 
must be continued. For practical purposes the search could 
start at any node which was not yet included in the search 
and which is identified by a "O" in LC. 

The newly found loops are recorded only if they 
are different from the ones already found. This condition 
can be detected easily, using logical functions of the pro- 
gramming language, by comparing the new jth entry row for 
matrix L with the other rows of L already set. 

The multiple starting of a search is necessary 
if the search was started in a part of the graph which is 
connected to the rest of the graph only by branches coming 


from the other part, as illustrated below in Fig. 5. 
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Bagi 2 


In the graph given (Fig. 5) a start at nodes 5, 
6, 7, or 8 could not lead to the detection of the loops 1-3-4 
and 3-4. On the other hand a start at node 2,1,3, or 4 would 
lead to the detection of all loops in one search. 

Here, how fast the search will be completed 
depends on the organization of the input information by the 


user. 


Ce THE PATH SEARCH 
ij tnibiabization 
For the purpose of finding all paths of the graph 
leading from a specified node to another specified node, a 
reorganized input list is formed. This list is called LIST2. 
As a first step a multiple node control vector 


LIST2(4) is formed, exactly analogous to LIST1(4). 
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2. - Pabhiseanon 

The path search follows the same logical sequence of 
events as the loop search, except that the path Search Stage. 
not at the element at the top of the list but with LIST2(l,a), 
using LIST2(3,a) as the first search key. The node "a" is 
the specified end node, to which a path is sought. 

The result of the search using the search key is an 
index k at which 

LIST2 (3 7a) =" DIsit2es ) 
Now the search key is changed: 

search key = LIST2(3,k) 
and the search continues as described before (B.2). This is 
the basic description of the path search. The difference , 
compared to the loop search, is the fact that here the pre- 
decessors of the end node are found and that here all connec- 
tions are avoided which would result in loops or which would 
not end at the desired starting node. 

a. Detection Of a Path 

A path is detected when it is found that the 

search key is set to an element LIST2(3,1), where the node 
"i" is the specified starting node. The path so detected 


a" to node "i", or, in the Notatron used 


leads from node 
here to facilitate computer programming, from LIST2(l,;a) "Ge 
LiST2.(\35)e 

During the search a temporary control vector is 
updated with "1" entries as described above (2.b) for the 


loop detection. Here the element LIST2(1,a) and all the 
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elements LIST2(3,n), which are used as search keys are marked 
with a "1" in the related place in the control vector. If 
for the element LIST2(3,i) a "1" is already marked in the 
control vector, a loop is detected. This is an unwanted 
Situation. The loop must be left by going back to the last 
multiple node. Then one of the other alternate paths going 
back from that node must be followed. 

How is it possible to keep track of the multiple 
nodes passed? It must be remembered, that the multiple nodes 
are marked in LIST2(4). (See definition for the analogous 
vector LIST1(4)). So every time a new search key is selected, 
LIST2 (4) will be checked, and a "1" in the position k, which 
is the "found" - index, indicates that the node LIST2(1,k) is 
a multiple node. When passing these nodes during the search, 
the element LIST2(4,k) is set to "0" so that it cannot be 
found again. 

The control vector BC (for Branch Control) is 
used to store all multiple nodes passed during the search, 
probably best in the form of the index, where the multiple 
node can be found in LIST2(1). It is then possible to 
identify the node to which one must go back to resume the 
search, in case a loop is detected, or the path detected ends 
at a node not specified as the end node. 

The latter case is possible and is identified 
by the fact, that using the present search key, (the last one 
assigned) no element LIST2(1,i) can be found in column LIST2 (1) 


which is equal to the search key. 
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b... Recording ©Ofsthe wath 

To record the actual sequence of nodes in the 
path, backwards from the end to the start, it is necessary to 
store in a temporary list all the elements LIST2(3,n) used as 
search keys during the search and the initial element LIST2(l1, 
a). This list contains only those nodes which actually form 
the path. Intermediate restarts of the search are not 
recorded. 

After the search for the jt path is completed, 
the jth row of the matrix P is filled with the appropriate 
entries. This can be most practically accomplished by 
defining the temporary list mentioned above in the form of a 


h 


row of P. Then it is only necessary to set the ;¢ row of P 


equal to this temporary list at the end of the search. 


De EVALUATION OF TRANSFER FUNCTIONS 
To evaluate the transmission gain from an input node "i" 
to an output node "o" (here: a node, which must have at least 
one incoming branch) the Mason Gain formula is used: 
LA, P 
[= k 
where P, is the gain of the kth girect path from node "i" to 


node "o 
Ay eethe c@iPaAcCcor OF che kth path, formed from 
A by striking out all” terms containing loops whe 


are touched by the kth patil. 
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-. evaluartion or the,Graphy Determinant? A” (DELTA) 

If the graph would contain only loops which do not 
touch, then the graph determinant would be evaluated as 

DEL tee Beeb) (l-lp). ~~. (1—lise) 
where i = nl, the number of loops in the graph, and L5 is the 
loop gain, the product of the gains of the branches which form 
the loop L.. 

If touching loops are present, then the products of 
their loop gains are eliminated from the expression for DELTA. 

a. Evaluation Of the Matrix LT 

The matrix LT contains information about which 
loops touch each other. 

The matrix LT is evaluated as: 

LT = L @ Lt L? is the transposed matrix L, 

@ is the operation defined in 

“eg oP 
Proof. Two loops touch each other when they have at least 
one node in common. By definition of the matrix L and the 
operation @ an element LT(k,j) = 1 if at least one of the 
logical expressions aj, & Dye = 1. There are k of these 
expressions, for every node in the graph. 

By definition of the operation " & " an expression 
alee & DK 5 = l only if a;, = 1 and also ae = J. But this is 
the case only if the node k is in loop i (for asy) and also in 
Hoop j (£65 by) Gensequently. DL1 ,jle= 1 only ui one of the 


graph is in loop i and loop j simultaneously. This completes 


the proof. 
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It should be noted that the matrix LT is 
symmetric and has always all elements on the main diagonal 
equal to l. As a consequence of this all the elements above 
or below the main diagonal and the elements on the main 
diagonal could be discarded. The computer program, Mason 
which is used to evaluate DELTA (Program 2) uses this fact. 


For convenience the example of Fig. 3 is given: 





O Oe 
Il 
1 i 
O 1 gill 





Now LT gives the result: 


L5 and L, touch, 
L, and Ly touche 
The number of touching loops in the graph (of yorden.2 jis 


equal to the sum of the elements of LT divided by two. 
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b. Evaluation of DELTA By a FORMAC Program 

A computer program, uSing the language FORMAC, is 
@eed to evaluate the expression for DELTA. (Program 2) 

The program is only meant as a demonstration of 
how this problem of evaluating DELTA can be solved using an 
algebraic manipulation language. 

After each multiplication by (1-L, ) the expanded 
result is checked for the presence of products of two touch- 
ing loops, and if such a product is found, it is replaced by 
zero. It is not necessary to check for products of more than 
two touching loops, because larger products containing the 
smaller one will go to zero automatically, after the smaller 
product is set to zero. 

To find the sets of two touching loops, the upper 
triangular part of LT, excluding also the main diagonal, is 
read row by row. The appearance of a "1" during that process 
Signals a touching of the respective loops, indicated by the 
row- and column index of the "1" found. This information is 
stored in a table in the form of an algebraic expression 
consisting of the product of the two loops involved. From 
this table the expressions are taken when, after the multi- 
plications mentioned above, a check for their presence in the 
accumulated product is made. 


2. Evaluation Of the kt® cofactor A, (DELTK) 





Having evaluated DELTA, the graph determinant, the 


kth cofactor is evaluated in a similar way, using the same 


3 


computer program. Instead of the matrix LT now the path- 
loop touching matrix PL is used as input. 

The matrix PL is evaluated as: 

PL = L @ ptr 
where P! is the transposed path matrix P. The resultant mat- 
rix PL shows the path - loop connections. If the matrix is 
read row-wise, then for each row i and column j it is indicated 
by a "1", if loop i touches path j. The proof follows the 


proof for the case of the touching loop matrix LT. 


For the example given for the loop case: (Fig. 3) 
J 
0 
PL = 
it 
1 





which indicates, that only loop 2 is not touched by path 2 
OL Daewoo 

Now two possible ways of solution can be followed. 
One way 1S to analyze the matrix PL in a way Similar to the 
one used for the matrix LT before. For each column of PL a 
specific cofactor DELTK, is evaluated. For a "1" in row i, 
loop L(i) is eliminated from the previously evaluated expres- 
sion for DELTA. Now the simple replacement of L(1) by zero 
results in elimination of all terms in DELTA which contain 
that loop L(i). After processing in that way all columns of 
PL, all cofactors are evaluated. Now DELTA and the evaluated 
cofactors DELTK, for the k paths present can be used in the 


final evaluation. 


Bh 


The other possible solution is first to use DELTA in 
the final evaluation by Substituting for every L; in DELTA 
the actual algebraic expression for L; and then evaluate the 
cofactors new from the Original expression involving the 
accumulated products of (I-L,), Similar to the evaluation of 
DELTA. 

3. Evaluation Of the kth Path Gain (PK) 

The algebraic evaluation of the path gain is best 
done immediately after a path is found (see II.C.2.b). One 
method of solution would be to store during the search for 
the path also the index k of the different Search keys. This 
index sequence could then be used in the algebraw@=part of 
the program to find those branch gains in the input table 
LIST! which must be multiplied to form the path gain. 

Two basic approaches are possible: 

~ the symbols stored in LIST1 (2) represent 
larger and more complex algebraic expressions. 
Then they constitute names or, to use FORMAC 
vocabulary, they are FORMAC variables. In this 
case there must be an assignment of the alge- 
braic expression to that variable somewhere in 
the program. (See example below) 

- the symbols stored in LIST1(2) are the actual 
algebraic expressions which constitute the 
branch gains. In FORMAC vocabulary they are 


called FORMAC expressions. In this case there 
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pale 16 


must be an evaluation which calculated the path 


gain directly. (See example below) 
Example. 
Case l. 

LIST1 (1) ErsemMr ) Eto Wl (3) 

1 b 2 

z a 1 

2 c 3 

2 k 6 

2 1 7 

3 e 4 

4 d 3 

4 g 5 

5 h 6 

6 a: 7: 

7 si 4 

7 g 6 


Note: the elements in LIST1(3) must be character type. 

In PL/1 this is indicated by writing for the character 
b the expression 'b'. (See PL/1 reference manual) 

A FORMAC statement like the following would produce 
the path gain in terms of the letters in Moers): 

bole) = "LIstL(2 1)" * “LeueeC2 3 he 

meee LIST) (2a) : 

The following statement would produce this output: 

PRINT Out (PK(1)); 

PIG) = Be Bn... 

Because generally it would not be known explicitly which 
symbols out of LIST1(2) form the gain PK, the following 
assign procedure would be used: 

At the time a more specific evaluation in terms of 


the actual algebraic expressions has to be produced, an input 


2) 


statement 1s executed. This statement would cause the read-in 
of the actual algebraic expressions. Then a replacement 
would take place, in which the letters a,b,... are replaced 
in the gain expression PK({1l) = BCEG... F by the algebraic 
expressions. Then a further expansion of the gain expression 
could take place. 

On data cards: 

'K7' §'3 * AS 2 Ss ee 

Then the statements: 

DECLARE EXPR CHARACTER(80) VARYING; 


to define EXPR a character variable, 
to define an index. 


K = 0; 
A: GET LIST(EXPR); K = K+l; 
ON ENDFILE(SYSIN) GO TO END; 
LET (PK(1) = REPLACE (PK(1), PLIST (2,K)" fe’ EXPR oe 
GO TO A; 
END: LET(PK(1) = EXPAND (PK(1)) ); 
PRINT OUT (PK(1)); 
The printout would read: 
PK(1) = K7 3 A* / 8s 
In this case the further expansion has no effect, 


because there are no expansions to be performed. 


Now consider case 2: 


LIST1 (1) List t2) LIST1 (3) 

as K7 as 

before 3*A** 2 before 
7S 
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In this case there are no data cards. The statements: 
Note: NB = number o- branches 
in a path 


LET (PK(1) = 1); 
DO K = 1 TO NB; 
LE rtp) = PK (ls YU LISTL(2,INX(K).)" ); 
END; 
PRINT OUT(PK(1)); 
Note: INX is an array of indices 


which indicate the branch 
to be included in the 
preduet. 
In the example given the array INX would contain for 
the first path PK(1l) the following indices: 
INX (1) 
INX (2) 
INX (3) 
INX (4) 


INX (5) 
INX (6) 


reo O NWF 


and NB = 6. 

A similar evaluation has to take place for the loop 
gains, the values of DELTA and DELTK. 

4. Final Evaluation Of the Transmission Gain 

The transmission gain G can finally be evaluated 
uSing a FORMAC program. Depending on the complexity of the 
problem first G could be evaluated in terms of the branch gains, 
where the branch gains are represented by their variable names, 


and by the loop gains. In this case G would be an algebraic 
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expression such as 


G,. <uamelalesl MUTI) gestenas) cele Co 
(ar (Ds) 


That is, DELTA and DELTK(1i) would be in terms of the 
loop gains; L(1i) and PK(i) would be in terms of the branch 
Gains apormc. .. 

Then after G is evaluated as above, a substitution 
could be made, where, again depending on the complexity 
involved in possible applications, several methods are possible. 
In this process the FORMAC functions REPLACE, CHAIN, EVAL and 
EXPAND[1L] would probably have to be used. Generally a com- 
promise has to be made between using more time for execution 


and less core space and vice versa. 
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III. DIRECT EVALUATION OF 
DETERMINANTS USING FORMAC 


In this chapter the direct evaluation of determinants 
using FORMAC programs is shown. By "direct" evaluation is 
meant a process of evaluation where the properties of deter- 
minants are used rather than graph theory as used in chapter 
II for the evaluation of DELTA (which is a determinant). 

Ae. A FORMAC PROGRAM USING THE BASIC DEFINITION OF THE 

DETERMINANT 

In numerical mathematics a determinant is normally eval- 
uated by normalizing the rows or columns and using elementary 
row- or column operations to eliminate all elements off the 
main diagonal. Then the value of the determinant is the 
product of the elements on the main diagonal. If this process 
would be used in the evaluation of determinants with algebraic 
terms as elements, the generation of rational functions or 
otherwise more complex algebraic expressions would create 
additional complexity. The handling of rational polynomial 
expressions, when the polynomials are in symbolic form, is very 
time consuming. 

The program DTERM uses the basic definition of the 
determinant: 

If Ais a square matrix of order n, then the determinant 
of A is DET(A): 


ce 
DET(A) = £(-1) ay a. ares | 


he 
J i ipa Jn 
where the second subscripts j assume all possible arrangements 


in which each column is represented exactly once in each term 
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of the sum, and the exponent, t, is the number of inter- 
changes necessary to bring the second subscripts into natural 
order (that is fel a2 a, 

In other words, all possible permutations of the numbers 
1,2,3, ..., n form the complete set of second subscripts 
mentioned above. 

If the number t is even, the permutation is called "even", 
otherwise "odd". 

The permutation of a set of numbers soon exceeds the 
time available for execution of any program. However, for 
n= 2 to n = 9 one may find it possible to use this method. 
The algorithm used in DTERM to find all the permutations of 
the numbers 1,2,3, ...-, n, for different n has the following 


execution time: 


n= 4: 0.05 sec 
n= 5: 0.09 sec 
n = 6: 0-31 sec 
n= 7: [ago e sec 
no = 8 “= ols. Zesee 


no = 0 eel oa See 

In the example program DTERM the permutation of column 
indices is only used to demonstrate the algebraic part of the 
program. Normally this algebraic part would be used with 
another method for finding the elements of the matrix A which 
have to be multiplied to form a term in the expansion of the 
determinant. For very large matrices, (for example, a 3/ x 
37 matrix) the graph analysis as described before would 
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identify those elements. Also the method used by Lipson [4] 
could possibly be used to identify those elements. How much 
computational effort each method requires in terms of time and 
Space has to be determined. 


1. General Description Of the Evaluation Of the 


Determinant By the FORMAC Program DETERM 
a. The Non-Algebraic Part Of DETERM 


The first step after the initialization is the 
evaluation of a permutation of the numbers 1,2 ..., n. The 
program used for that purpose is the subroutine PERM’. The 
program PERM produces upon each call another distinct permuta- 
tion of the n ! possible permutations in such a way that the 
permutations it generates are alternately odd and even, a 
property the author of algorithm 115 PERM described as "prob- 
ably useless". However, this property enables us to determine 
the sign of the permutation. 

Having found a specific permutation it is used as 
the set of indices jl, 32, ..., jn mentioned in the definition 


of the determinant. At the same time the sign of the product 


a a | is known. 
a 
dh Z n 
The elements a. in the product are then tested 
J 
for being zero. If a zero element is found, a new permutation 


is produced and the process starts again, because further 


algebraic manipulation would only produce a zero term. 


1fhe authors conversion of the algorithm 115 PERM, H. F. 
Trotter, Comm. ACM, Aug. 1962 from ALGOL into PLl. 
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In case the elements are all non-zero the product 
is algebraically evaluated. This is a complicated step, 
which will be described in the next section. 

Then the next non-zero term (a term is one of 


the products as. 


ees aaa 1S evaluated in the same manner. 


Finally all coefficients of equal powers of the 
polynomial variable in the different terms of the determinant 
are collected. 

b. The Algebraic Part of DETERM 

Most of the algebraic manipulation statements 
are. detectable bygthe word “rer Qi | ee Se ); and the 
embracing brackets. The semicolon is a feature of PL/I) ae 
Signals the end of a statement. 

Consider now two polynomials in the variable S: 

Pl = (A + B*)2 § + (C + D(E + F)3 s2 

P2 = (A+C)? § +B S3 + F(G-H) 85 
If the product P1*P2 is formed and all the terms are completely 
expanded, a new polynomial is formed which has 90 terms before 
Simplification. 

If the same product is formed without expanding 
the coefficients of the variable S the resultant polynomial 
has only 6 terms. Each of the coefficients of the variable § 


is in unexpanded form. 


Example. 
Ple=» (a +eb)esm+ qemqeeqpeeg? 
P2k= (ees age 
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The expansion used in DTERM: 
Pl P2.2%(a + bye + fF) s4 #%(a tb) (g +h) s 
+ (c + d)(e + £) s? + (c + a)(g +h) s? 
Then this result is simplified to 
Pl P2 = (a + b) (g + h) s + 
Glas be (e + £) + (c + d) (g + h) ) Ss? + 
(c + d)(e + £) s? 
or = As +Bs2+4+C s3 
where the coefficients A, B, C ... are stored in this form on 
disks. After all terms of the determinant are expanded in 
this manner, coefficients of equal powers of s are read back 
into main core and are expanded fully and added. They are 
then printed and the space used by them is subsequently made 
available to other uses. 

The latter method is used in the first part of 
the algebraic manipulation to save space. A main limitation 
of FORMAC (probably of any algebraic manipulation language) is 
the enormous space requirement for certain expansions; there- 
fore to make algebraic manipulation of large expressions 
possible at all, it must be done at the expense of time. It. 
has been shown [6] that running FORMAC programs can be made 
much faster after some refinement; often the omission of the 
use of one FORMAC function showed a considerable increase in 
execution time. The reason is that each FORMAC function use 


constitutes really the use of a subroutine. 
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The main emphasis in DETERM was placed on the 
possibility of evaluating very large determinants, the 
refinement was secondary. 

After all n factors, where each is considered 
being a polynomial, have been multiplied, the resultant 
product, which is called "PD& inethesprogram, is furthes 
processed. 

Each coefficient in PD of the variable S and the 
constant coefficient is identified and stored separately on 
disk file. Then the evaluation of a new term PD starts. 

After all terms PD, which are the products 
mentioned in the definition of the matrix, have been processed 
in the manner described above, the coefficients stored on disk 
are now retrieved and added up, always adding coefficients of 
the same power of the variable S originating from the different 
terms of the determinant. The coefficients are designated as 
follows: 

COEF (1).=.¢C(1,1)+ C (233) eee 1) 
where i is the power of S which is under consideration and the 
numbers 1,2, ...k are the numbers of the k non-zero terms in 
the determinant. 

The coefficients C,; are the ones which are 
stored on disk. The final coefficient COEF; is the coefficient 


of the variable power si in the expanded determinant. 
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Example for using DTERM? 


System equations: 


e. =k . 6. 
1 pi 2 
e =k 8 

O pO oO 
e =e, - =n 
st . Me a 
e, = In! a 
as = = - e, 


E. (s) sta 


Eas). Con Sc aaa 


The state equations are formed with the following variables as 





state variables: 


xX, = 8 A 

2S = oe 

X3 = i. 

Xi = dummy variable 

Bic) Xm Ee i sta 





“Example taken from notes course 3411 Naval Postgraduate 
School. 
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which gives: 


The equations in the S - domain: 


CaaS 


Zr {104 


(3) — Rays! a (ko h 1/p.) Xo + (r-) X3 
+ (Ls) X3 + (k_ (b - a)) X 
(4) (s) X 


The equations in matrix form: 


AA = 
where A is 
X is 


B is 


The elements 


Sawai 


ee 


Ag? = 


Ay, 


2h eS 


i dee 


+ (k ) X 
po 


4 il 


B 


a matrix of polynomial coefficients 
a matrix of variables 


a matrix of constants 


2 + Jm) s Xo + (f£m(p, ag) + 


+ (a + b) X 


of A as used as input to DTERM: 


(jo (py o) 4 + Jm) s + (fm(p, ois + fm) 


-— (P17 Po Ky 1,) 
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k_k_. 
a pi 


i 


eal 7 Oe sey 


ie 


A, = (Ls + r.) 
A = (k 

34 ( - 
Ang = (s t+ at b) 


The program DTERM with the matrix A as input and the output 
1s contained in section COMPUTER PROGRAMS. 
Figure 8 shows the physical system, Fig. 9 shows 


the Mason flow graph. 
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49 


2. Discussion Of the Flow Graph For DETERM 
The flow graph for DETERM is discussed here, to some 
limited extent. 
For clarification the definition of the determinant 
is now again given in terms of the symbols used in the program. 
DET(A) = PD, + PDot+ . . . # PD. 
where 
PD, = C(3,i) si i = 0,1,2 
and finally 
DET(A) = COEF(0) + COEF(1) S + COEF(2) S*# +... 
+ COEF (MxPow) sMAPOW 
where 
COEF(i) = €(1,2) + C( 23g .:. weGeny.) 
a. Data Design For DTERM 
Consider now the computer program printout. 
(Program 1). All statements beginning with "DCL" are declara- 
tions, they define those data elements used in the program 
which are defined explicitly. Some of these definitions will 
now be discussed. 
IP is an integer array. It contains after Gach 
call of PERM a distinct permutation of the numbers 1,2,...n. 
These numbers are used as column indices in designating elements 
of the input matrix A. The array Albis a matrix of characters, 
Each element can have a varying length of maximal 75 characters. 
PDD is a character type variable. It contains at 


certain moments an element of the matrix A.N1 is the order of 


the input matrix A. 
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General flowgraph of DTERM 


PART J 
PART 2 


find ao permutation J 
of column indices 











Is 
the ith 
product =O 










Evaluate the product PD. 
The coefficients C; are not 


expanded. 













PART ¥ 
Extract the coefficients Cj 
from PD;. Store on disk us- 
ing SAVE. 

PRRT S 






increase 
index 






YES 


Add all coefficients of equal 
powers of the variable s. PART 7 


This results in C (i). 
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SPECIALIZED FLOWCHART OF DTERM 


DTERM 
Definition of data elements PRRT L 


Initialization and input of matrix A PART 2 
C [index | 


Determination of sign of the (F) 
permutation i. 


Determination of permutation |. 


Increase index IN IN Counts the number of 


nonzero products in the 


Initialize for product accu- index M | 


mulation of PD. 











PART 3 


Increase 
Index | 


) 
€ 


Use the jth argument of new 
ndex element as multiplier. PART ¥ 
A) 


@ 
Increase index 





Multiply the Ith argument 
of the previous result PS 
by multiplier. 







Result of this step is assigned 


another name to moke it the PART ¥ 


‘next ‘previous’ result. PD=P5 





‘> NO 


YES 


Find highest power of S in all 


products ein” MAXPOW. 





PRRT S 


Find all coefficients C(i) of s' 
PD . Store on disk using SAVE. 










PART 6 





nonzero products 
of DET (A) found 


NO 


( ) YES 


Take coefficient C., (i) of s'in git 
eocedl |i Per? 






Mis arunning index 


is arr (i from 7 to IW. 






defined 


YES 
Hy 
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Determine the exponent of the 


vorioble s for thot result. 













Is 
exponent = O 
(result is 


YES Add result of 
multiplication 





















const. ) to previous 
PART ¥ result P5. 
Find coeffient of 
same power in pre- Store PS on 
vious result P5. disk using 
SAVE. 
Is 
this 
sc coefficient 
present 


in P5 


Add result of multiplication to 
previous result P5. 









Build an argument 
list for REPLACE: Store P5 on disk 
function. using SAVE. 


Replace the old coefficient of PART ¥ 


S we with the sum of the old 


one + result of multiplication. 






All 
orgumenis 
are multiplied 


YES NO 
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& 














ls 
there o 


previous sum 
? 


NO 


COEF (i)=C,, (i) 






COEF (i) = previous sum 
+ new Cy (i) 


‘ PRRT 7 


Print COEF (i) 


NO 
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b. The Logic Sequence (See flow chart, Figq-ey 

The initialization includes some statements 
which are self explanatory. Some will be discussed now. N2 
is used in a loop to fill the array IP initially with the 
sequence of numbers 1,2, ,,,n in natural order, an input 
requirement for the subroutine PERM. Also the binary variable 
"FIRST tseusedefor comtrol of PERM. The following Woere 
uSing the index I=l to N! is the main loop of the program. In 
that loop all n ! possible products of the determinant are 
investigated. Then the sign of the permutation is determined. 

The loop M=l1 to N=l includes the zero check for 
the ith product of n !_ products. 

The next loop M=l1 to Nl is reached only if all 
factors in a product are non-zero. Then the algebraic manipu- 
lation begins. The procedure was described in 2. befom. 
The action of the various FORMAC functions (NARGS, ARG, LOWPOW, 
ATOMIZE, etc.) is described in the FORMAC reference manual [6]. 
Bs EVALUATION OF A DETERMINANT USING BARREISS FRACTION-FREE 

ALGORITHM 

Another method of evaluating a determinant is using an 
algorithm given by Lipson. This algorithm is based onvageva. 
step variant by Barreiss of a fraction-free elimination algo- 
rithm, which in turn is a modification, attributed to Jordan, 
of Gaussian elimination [4]. 

Lipson uses the algorithm to evaluate a system of linear 


algebraic equations by transforming a matrix into triangular 
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form and following back-substitution. The whole process is 
fFraction-free. 

The linear equations processed have symbolic coefficients 
in the form of multivariate polynomials. 

Lipson has also programmed the algorithm in the language 
"SCOPE FORMAX". [4]. 

1. The Fraction-Free Gaussian Elimination 

Consider the n x (n+m) augmented matrix of a system 

of linear equations Ax = B: 


Aj 4 (0) = (ais (0)) = (A|B) 


Then the modified algorithm is given 


afk-1) g(k-1) - (kel) (kel) 


.,, s kk ij j ik 
1) q \K-2) 
k-1,k-1 
where 
Ker = “lecyence, nol 
ee es sp TL 
foie, 6 > NFM 


and with a‘) as given before, and where 


al) 
ang = 1 
lle i =1,-.-,k; 4 <= 1,-+.,nt+m 
ame = 
J 0 a ea Sie: are 


Also, given here without proof (for the proof see 


Ref. 4) 


5) /| 


The determinant of the matrix A 


(n=) 
nn 


DET(A) =a 


The formula given above was programmed using 
normal FORMAC, also ignoring for purpose of clarification any 
pivot element check for zero. 

The algorithm alone does not give a fraction-free 
algebraic representation. Some processing is necessary to 
bring dividend and divisor into a form such that cancellation 
takes place. 

A very useful method may be to use this algorithm 
to compute the determinant of a matrix with single letters as 
elements, followed by a substitution process, where more 


complicated expressions are substituted for the letters. 
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IV. APPLICATION OF FORMAC 
IN OTHER ALGEBRAIC PROBLEMS OF CONTROL THEORY 


In this chapter some other applications of the algebraic 
manipulation language FORMAC in control theory are proposed. 
Some examples will be given, the programs being more of demon- 
stration programs than operational programs. 

The reason for this is mainly the lack of time which is 
needed to prepare an optimal program, optimal with respect to 


execution time and space requirement. 


A. ROOT LOCUS THEORY 

The impulse response of a linear time-invariant system is 
determined by a sum of exponential time functions. The expo- 
nents of these terms are the roots of the characteristic 
equation of the system. 

A system is stable if and only if none of the n roots of 
the characteristic equation has positive real parts. 

Consider the polynomial in the Laplace-variable S which 


represents a characteristic equation: 


The coefficients a, can be complex. Some of the coeffi- 
cients can include a design parameter, Kk. 


The polynomial F(s) can be partitioned into two parts as 


F(s) + K B(s) (1s) 
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The complex variabWe S 1s replaced as™s — 0 + J 0. ie 


require the real and imaginary parts of F( o , w ) to be zero 


independently: 
Re A(a,w) + K Re Bl(o,w) = 0 (2a) 
Im A (GRP Kk Ina (ou) —= 0 (2b) 


Using the set of simultaneous equations (2a) and (2b) the 
parameter K can be eliminated: 

Re B * Im A + Re A * Im B = 0 (3) 

This equation is the equation of the root locus; that is, 
every root of the characteristic equation F(s) = 0 muSt Sacuee 
the equation (3). 

In the computer program (Program 3) this equation is 
called RLOC = f(0,w). RLOC contains in every term the 
variable w, because the real axis iS always a locus for the 
roots. This factor of w can be removed from RLOC. 

If now a value for one of the variables o or w 1s 
assumed RLOC will give those points in the s-plane which 
satisfy equation (3). These points are on the root locus. 

Il. A FPORMAC Program FOr vEnewRooEmuocuS 

As an example for the application of FORMAC a program 
is given (Program 3). 

The program reads in a polynomial CHEQ and the para- 
meter PARAM in algebraic form. Then the polynomial is expand- 
ed and then, following the development as described before, 
the equation (3) of the root locus is derived in algebraic 


fOr Ms 
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Ac this@poeint values fer%Gror o°in the range “or 
interest can be chosen and the equation for the root locus 
canmbe evaluated, givangevalues fer wero respectively. 

Wormally the expression in the remaining unknown 
will not be linear. In that case the expression is treated 
aS a polynomial, the roots of which can be found uSing a 
Boot-—fanding groutme. 

Finally the results can be plotted, giving the 
complete root locus in the range of interest. 

To evaluate the values of the parameter at points 
along the root locus it is necessary to substitute back the 
values of o and w for those points into equation (2). 

2. Discussion 

The conventional root-locus procedure [9] is basically 

a trial-and-error method: 


For a range of parameter values Ky 7Ko,---, K POX 


ne 
which one suspects an acceptable area of root locations in 

the s-plane, the roots of the characteristic equation are 
found by a root-finding method, probably a computer-programmed 
routine, uSing a standard method. 

The the root locus is inspected in the area of 
interest for the values of o and w and the corresponding para- 
meter values are selected. 

The problem is that values of a quantity which are 
sought have to be assumed beforehand. 

The method proposed here [10] avoids this sequence. 
The equation of the root locus is evaluated in the range of 
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interest for the variables o and w. Then by back-substitution 
of o and w into the given characteristic eGureton the para- 
meter values are calculated. Another advantage is the fact 
that the degree of the root locus equation (3) is generally of 
lesser degree than the characteristic equation. This face 
results from possible cancellations of terms of equal powers 
in the process of parameter elimination. In example 1 the 
degree of the characteristic equation is 4 while the equation 
of the root locus is of degree 3 ino and of degree 2 in w. 

Bs FREQUENCY FUNCTIONS 

Consider a transfer function P(s) = Y(s)+X(s). P(s) 
might contain some parameter in symbolic form. It is possible 
to use a FORMAC program to evaluate the absolute value of P(s) 
and the angle of P(s) as expressions containing the parameter 
in symbolic form. For complicated transfer functions 
might be easier to evaluate the complex algebra only once and 
then vary the parameter to obtain values for magnitude and 
angle of the transfer function. 

For the latter purpose it is mentioned that in FORMAC 
transcendental functions with symbolic arguments can be used. 
Cs LAPLACE TRANSFORMS AND INVERSION OF LAPLACE TRANSFORMS 

Consider a differential equation in the time domain. If 
such an equation is given as input to a FORMAC program it is 
possible to perform a Laplace transform on the equation. The 
basic problem to be solved would be the replacement of one 


@eematon by afidther one. “9S boehnmoides oretherditterential 
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equation would be read in as two separate FORMAC variables, a 
transfer function could be formed. 

Another definite possibility is the analysis of an alge- 
braic expression in the variable s and subsequent recognition 
of certain algebraic patterns, which can then be replaced by 
an expression in the time domain. This method would be 
Similar to the "manual" approach, because one often solves 
these problems in the same way: recognizing the transform 
pair from experience (or tables). 

De SENSITIVITY ANALYSIS 

A very powerful feature of FORMAC is its ability to form 
partial derivatives of user specified order. Also the mathe- 
mMatical functions such as sine, cosine, exponential, logarithm 
can be used with symbolic arguments. 

Consider a mathematical model T(k), which might be a 
transfer function or frequency response function, of a linear 
time-invariant system. 

The sensitivity of T(k), k being a variable, with respect 


to k is defined: [9] 


Aikido es) GUTKk) ok 
k “dink sae) | «TK) 


iiewoctorelvrety or the phase angle of T(k), 9%, with 


respect to parameter k is defined: 


3 fT _ dep k 
k dk mp 


A FORMAC program which would find the derivatives involved in 


the sensitivity analysis could be time saving, if symbolic 
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arguments were involved and one would like to vary the para- 
meter. A program for the finding of the derivatives would 


only be used once. 
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V. CONCLUSION 


A method of solution of a pattern recognition problem in 


the analysis of a Mason flow graph has been shown. It also 


has been shown that computer programs, uSing an algebraic 


manipulation language FORMAC can be used in the subsequent 


evaluation of the results of the pattern recognition problem. 


It 1s possible to evaluate in symbolic form certain system 


functions, for example the transfer function between specified 


nodes and the characteristic equation of large systems. 


A FORMAC program can also be used to advantage in a 


modified root locus method. 

It is felt that the application 
lation language (FORMAC) which has a 
as compared to other more restricted 
handle only polynomials, is of great 
algebraic problems in control theory 


Of particular importance is the 


of an algebraic manipu- 
wide range of capabilities 
languages which may 
value in the solution of 
and in network theory. 


fact that FORMAC isa 


Superset of PL/1. All the capabilities of PL 1 are usable 


together with FORMAC. 
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COMPUTER PROGRAMS 

Program DTERM, a FORMAC program 

Purpose: To evaluate a determinant of a square matrix 
of order up to 9. The elements of the matrix 
are univariate polynomials with symbolic 
coefficients. Given is one example. 

Program FINDLOOP, a PL/1 program 

Purpose: To evaluate a control vector LIST1(4) which is 
needed for the Mason flow graph analysis. Given 
is one example. 

Program MASON 1, a FORMAC program 

Purpose: To find sets of touching loops from given topo- 
logical matrices describing a Mason flow graph 
and to evaluate the graph determinant. Given is 
one example. 

Program MASON 2, a FORMAC program 

Purpose: To evaluate the cofactors in the Mason gain for- 
mula. Given 1s one example. 

Program RTLOC, a FORMAC program 

Purpose: To find the equation of the root locus from a 
given characteristic equation in symbolic form. 


Given are three examples. 
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